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Tsallis' q index and Mori's q phase transitions at edge of chaos 
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We uncover the basis for the validity of the Tsallis statistics at the onset of chaos in logistic maps. 
The dynamics within the critical attractor is found to consist of an infinite family of Mori's q- 
phase transitions of rapidly decreasing strength, each associated to a discontinuity in Feigenbaum's 
trajectory scaling function a. The value of q at each transition corresponds to the same special 
value for the entropic index q, such that the resultant sets of g-Lyapunov coefficients are equal to 
the Tsallis rates of entropy evolution. 
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I. INTRODUCTION 

Searches for evidence of nonextensive , properties 
at the period-doubling onset of chaos in logistic maps - 
the Feigenbaum attractor - have at all times yielded affir- 
mative responses, from the initial numerical studies 0, 
to subsequent heuristic investigations 1, and the more 
recent rigorous results 0,0. However a critical analysis 
and a genuine understanding of the basis for the validity 
at this attractor of the nonextensive generalization 
of the Boltzmann-Gibbs (BG) statistical mechanics - 
here referred as q-statistics - is until now lacking. Here we 
clarify the circumstances under which the features of q- 
statistics are observed and, most importantly, we demon- 
strate that the mechanism by means of which the Tsal- 
lis entropic index q arises is provided by the occurrence 
of dynamical phase transitions of the kind described by 
the formalism of Mori and colleagues 0. These tran- 
sitions, similar to first order thermal phase transitions, 
are associated to trajectories that link different regions 
within a multifractal attractor. The onset of chaos is 
an incipiently chaotic attractor, with memory preserv- 
ing, nonmixing, phase space trajectories. Because many 
of its properties are familiar, and well understood since 
many years ago, it is of interest to explain how previous 
knowledge fits in with the new perspective. 

The Feigenbaum attractor is the classic one- 
dimensional critical attractor with universal properties 
in the renormalization group (RG) sense, i.e. shared 
by all unimodal (one hump) maps with the same degree 
of nonlinearity. The static or geometrical properties of 
this attractor are understood since long ago - [To| . 
and are represented, for example, by the generalized di- 
mensions D(q) or the spectrum of dimensions f(a) that 
characterize the multifractal set 0, E3- The dynam- 
ical properties that involve positions within the attrac- 
tor also display universality and, as we see below, these 
are conveniently given in terms of the discontinuities in 
Feigenbaum's trajectory scaling function a that measures 
the convergence of positions in the orbits of period 2™ as 
n — > oo 8]. Let us first recall that the Feigenbaum at- 
tractor has a vanishing ordinary Lyapunov coefficient Ai 



and that the sensitivity to initial conditions £t does not 
converge to any single- valued function and displays fluc- 
tuations that grow indefinitely [l^, 0, 0, 0. For 
initial positions at the attractor £ t develops a univer- 
sal self-similar temporal structure and its envelope grows 
with t as a power law 0, [3, [3, 0, 0- We are 
interested here in determining the detailed dependence 
of the aforementioned structure on both the initial posi- 
tion Xo and the observation time t as this dependence is 
preserved by the infinitely lasting memory. Therefore we 
shall not consider the effect of averagin g wi th respect to 
xo and/or t, explored in other studies jl^ [l5| . as this 
would obscure the fine points of the dynamics. 

The central assertion of the q-statistics with regards 
to the dynamics of critical attractors is a sensitivity to 
initial conditions £t associated to the ^-exponential func- 
tional form, i.e. the 'q-deformed' exponential function 
exp g (x) = [1 — (q — Vjx] -1 '^ 1 ' . From such £ t a q- 
generalized Lyapunov coefficient X q can be determined 
just as Ai is read from an exponential £ t . The A g is pre- 
sumed to satisfy a g-generalized identity X q — K q 
where K q is an entropy production rate based on the 
Tsallis entropy S q , defined in terms of the q- logarithmic 
function ln 9 y = (y 1-9 — 1)/(1 — q), the inverse of exp q (x). 
Unlike Ai for (ergodic) chaotic attractors, the coefficient 
X q is dependent on the initial position xq and therefore 
X q constitutes a spectrum (and also K q ) that can be ex- 
amined by varying this position. 

The fixed values of the entropic index q are obtained 
from the universality class parameters to which the at- 
tractor belongs. For the simpler pitchfork and tangent 
bifurcations there is a single well-defined value for the in- 
dex q for each type of attractor as a single g-exponential 
describes the sensitivity [l9|. For multifractal critical 
attractors the situation is more complicated and there 
appear to be a multiplicity of indexes q but with pre- 
cise values given by the attractor scaling functions. As 
shown below, the sensitivity takes the form of a family 
of interweaved q-exponentials. The q-indexes appear in 
conjugate pairs, q and Q = 2 — q, as these correspond to 
switching starting and finishing trajectory positions. We 
show that q and Q are related to the occurrence of pairs 
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of dynamical 'g-phase' transitions that connect qualita- 
tively different regions of the attractor [l4|. These 
transitions are identified as the source of the special val- 
ues for the entropic index q. For the Feigenbaum attrac- 
tor an infinite family of such transitions take place but 
of rapidly decreasing strength. 

In the following section we recall the essential features 
of the statistical-mechanical formalism of Mori and col- 
leagues to study dynamical phase transitions in at- 
tractors of nonlinear maps and follow this by a summary 
of expressions of the g-statistics. Then, in subsequent 
sections we present known properties and develop oth- 
ers for the dynamics within the Feigenbaum attractor. 
Amongst these we derive the sensitivity £t in terms of 
the trajectory scaling function er, and use this to make 
contact with both Mori's and Tsallis' schemes. We dis- 
cuss our results. 



II. STATISTICAL MECHANICS FOR CRITICAL 
ATTRACTORS 

During the late 1980's Mori and coworkers developed 
a comprehensive thermodynamic formalism to character- 
ize drastic changes at bifurcations and at other singular 
phenomena in low dimensional maps 0. The formalism 
was further adapted to study critical attractors and was 
illustrated by considering the specific case of the period- 
doubling onset of chaos in the logistic map , 0] , . 
For critical attactors the scheme involves the evaluation 
of fluctuations of the generalized finite-time Lyapunov 
coefficient 



X(t,x ) 



t-i 

-^Vln 
\nt ^ 

i=0 



dfnoo( x i) 



dxi 



t > 1, 



(1) 



where /^(x) is here the logistic map, or its extension to 
non-linearity of order z > I, 



U(x) 



1 



•1 < x < 1 , < n < 2. (2) 



We denote by fiodz) the value of the control parameter 
H at the onset of chaos, with /!oo(2) = 1.40115.... Notice 
the replacement of the customary t by In t in Eq. , as 
the ordinary Lyapunov coefficient Ai vanishes for critical 
attractors, here at /i^, t — > oo. 

The density distribution for the values of A, at t ^S> 1, 
P{X,t), is written in the form Q 



P(X,t) = t-^p(o,t), 



(3) 



where ip(X) is a concave spectrum of the fluctuations of A 
with minimum ?/>(0) = and is obtained as the Legendre 
transform of the 'free energy' function <fi(q), defined as 



0(q) 



lim 



lnZ(t,q) 



t— >oo Int ' 
where Z(t,q) is the dynamic partition function 

-(q-l)A 



Z(t,q) = / d\ P(X,t) t 



(4) 



(5) 



The 'coarse-grained' function of generalized Lyapunov 
coefficients A(q) and the variance v(q) of P(X,t) are 
given, respectively, by 



A(q) EE *M and v(q) EE ^ 



dq 



dq 



(6) 



[T||. Notice the special weight t'^~^ x in Z(t,q) and 
in the quantities derived from it. The functions 0(q) and 
t/>(A) are the dynamic counterparts of the Renyi dimen- 
sions D(q) and the spectrum /(5) that characterize the 
geometric structure of the attractor . 

As with ordinary thermal 1st order phase transitions, 
a "g-phase" transition is indicated by a section of lin- 
ear slope m c = 1 — q in the spectrum (free energy) 
i/>(A), a discontinuity at q = q in the Lyapunov func- 
tion (order parameter) A(q), and a divergence at q in the 
variance (susceptibility) v(q). For the onset of chaos at 
fJ-<x>(z = 2) a single g-phase transition was numerically 
determined - [H| and found to occur at a value close 
to m c — — (1 — q) ~ —0.7. Arguments were provided 
- ^ijj for this value to be m c = — (1 — q) = — In 2/ In a = 
—0.7555..., where a = 2.50290... is one of the Feigen- 
baum's universal constants. Our analysis below shows 
that this initial result gives a broad picture of the dy- 
namics at the Feigenbaum attractor and that actually an 
infinite family of g-phase transitions of decreasing mag- 
nitude takes place at /ioo. 

Independently, Tsallis and colleagues proposed that 
for critical attractors the sensitivity to initial conditions 
£ t (defined as ^t(xo) = limAa: ^o(Axt/Axo) where Aa; 
is the initial separation of two orbits and Ax t that at 
time i), has the form 



£,t(x ) = exp q [Xg(x ) t], 



(7) 



that yields the customary exponential £ t with Ai when 
g — > 1. In Eq. J7J) g is the entropic index and the 
initial-position-dependent X q (xo) are the g-generalized 
Lyapunov coefficients. Tsallis and colleagues also sug- 
gested Q that the identity K\ = Ai [l^. where the rate 
of entropy production K i is given by 



Kit = Si(t) - Si(0), 4 large, 



and 



Si 



^2 Pi In Pi 



(8) 



(9) 



for an ensemble of trajectories with instantaneous distri- 
bution pi, would be generalized to K q = X q , where the 
g-generalized rate of entropy production K q is defined via 



K q t = S q (t) - S q (0), t large, 



(10) 



and where 



s q ^Y,p> ln op* 1= q _\ Pi ( n ) 
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is the Tsallis entropy. These properties have been cor- 
roborated to hold at ^oo(z) numerically 0, f° r sets 
of trajectories with x spread throughout — 1 < x < 1 
and analytically for specific classes of trajectories start- 
ing near xq = and observed at specific times of the form 
t = (2k + 1)2™ - 1, k = 0, 1, 2... and n = 1, 2, ... 0] - (H). 
We explain the rationale for these particular choices of 
xq and t in the following section where we examine the 
structure of trajectories inside the attractor. See Fig. 1. 



cally, the positions for the main subsequence k = 0, that 
constitutes the lower bound of the entire trajectory (see 
Fig. la), can be identified to be |af2«| — ^n,o = a - ™, 



where d n o = 



(2" 



'(0) 



is the 'n-th principal diameter' 

defined at the 2™-supercycle, the distance of the orbit po- 
sition nearest to x = [8j . The main subsequence can be 
expressed as 



exp 2 _„(-A 9 <) 



(12) 



III. DYNAMICS WITHIN THE FEIGENBAUM 
ATTRACTOR 

In Fig. la we have plotted (in logarithmic scales) the 
first few absolute values of iterated positions \xt\ of the 
orbit at (J,oo(z = 2) starting at xq — where the labels in- 
dicate iteration time t. Notice the structure of horizontal 
bands, and that in the top band lie half of the attractor 
positions (odd times), in the second band a quarter of 
the attractor positions, and so on. The top band is elimi- 
nated by functional composition of the original map, that 

(2) 

is by considering the orbit generated by the map (0) 
instead of /^(O). Successive bands are eliminated by 

considering the orbits of (0), j = 1,2.... The posi- 
tions of the top band (odd times) can be reproduced ap- 
proximately by the positions of the band below it (times 
of the form t — 2 + 4n, n = Q, 1, 2...) by multiplication 
by a factor equal to a, e.g. \x\\ ~ a|x2|. Likewise, 
the positions of the second band are reproduced by the 
positions of the third band under multiplication by a, 
e.g. X2 — a|x4|, and so on. In Fig. lb we show a 
logarithmic scale plot of 1— \xt\ that displays a band 
structure similar to that in Fig. la, in the top band 
lie again half of the attractor positions (even times) and 
below the other half (odd times) is distributed in the 
subsequent bands. This time the positions in one band 
are reproduced approximately from the positions of the 
band lying below it by multiplication by a factor a z , e.g. 
1 — |x3 1 ~ a z (l — \xs |). The properties amongst bands 
of iterate positions merely follow from repeated compo- 
sition and rescaling of the map and represent a graph- 
ical construction of the Feigenbaum RG transformation 
Rf(x) = af(f(x/a)), where a(z = 2) = 2.50290... 

Also, the trajectory with initial condition xq = maps 
out the Feigenbaum attractor in such a way that the 
absolute values of succeeding (time-shifted r = t + 1) po- 
sitions \x T \ form subsequences with a common power-law 
decay of the form t -1 / 1 ^ 9 with q = 1 — In 2/ lna(z), with 
q ~ 0.24449 when z — 2. See how positions fall along 
straight diagonal lines in Fig. la. That is, the entire 
attractor can be decomposed into position subsequences 
generated by the time subsequences r = (2k + 1)2™, each 
obtained by running over n — 0, 1, 2, ... for a fixed value 
of k = 0, 1, 2, ... Noticeably, the positions in these subse- 
quences can be obtained from those belonging to the 'su- 
perstable' periodic orbits of lengths 2™, i.e. the 2™-cycles 
that contain the point x = at JL n < /u c (0) 0. Specifi- 



with A g = (z — 1) In a/ In 2. Interestingly this analytical 
result for \xt \ can be seen to satisfy the dynamical fixed- 
point relation, h(t) = ah(h(t/a)) with a = 2 1 /( 1 ~ q \ See 
lU, @ for z = 2 and [HI for general z > 1. 

We now work out the relationship between the tra- 
jectory scaling function a and the sensitivity £j at /loo. 
To begin with we recall |S| the general definition of the 
diameters d n , m that measure the bifurcation forks that 
form the period-doubling cascade sequence. The d n . m 
in these orbits are defined as the distances of the ele- 
ments x m , Tn = 0, 1, 2, 2 n — 1, to their nearest neigh- 



bors f£ \x r . 



d = f (r 



h2"- 



>(0)-/J>>(0). 



(13) 



For large n, d n o/rfn+i ,o — —a(z); a(2) = a. Further, 
Feigenbaum [ll| constructed the auxiliary function 



cr„(m) 



d„ 



+ 1,771 



dn.7 



(14) 



to quantify the rate of change of the diameters and 
showed that in the limit n — > oo it has finite (jump) 
discontinuities at all rationals of the form m/2™ +1 . 
So. considering the variable y = m/2™ +1 one obtains 
[Tl| omitting the subindex n, cr(0) = — 1/a, but 
<r(0 + ) = 1/a 2 , and through the antiperiodic property 
a(y + 1/2) = -a(y), also a (1/2) = 1/a, but <r(l/2 + 
+ ) = — l/a z . Other discontinuities in a(y) appear at 
y = 1/4, 1/8, 3/8, etc. In most cases it is only necessary 
to consider the first few as their magnitude decreases 
rapidly. See, e.g. Fig. 31 in Ref. 0. The discontinuities 
of <r n (m) can be suitably obtained by first generating the 
superstable orbit 2°° at fi^ and then plotting the position 

differences \x t — x m \ — fj^ (0) — fj^ (0) for times of 

the form t = 2" — m,n= 0, 1, 2, in logarithmic scales. 
The distances that separate positions along the time sub- 
sequence correspond to the logarithm of the diameters 
d n ,m- See Figs, la and lb where the constant spacing 
of positions along the main diagonal provide the values 
for lnd nj o — — nlna and lnd ni i ~ — nlna z , respectively. 
The constant slope s m of the resulting time subsequence 
data is related to a n (m), i.e. a n (m) — 2 Sm . See Figs, 
la and lb where the slopes of the main diagonal subse- 
quences have the values — In a/ In 2 and —z In a/ In 2, re- 
spectively. From these two slopes the value of the largest 
jump discontinuity of a n (m) is conviniently determined. 
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FIG. 1: a) Absolute values of |a?t | vs t in logarithmic scales for the orbit with initial condition xo — at fi^, of the logistic map 
z = 2. The labels indicate iteration time t. b) Same as a) with \x t \ replaced by 1 — \xt\. 



A key factor in obtaining our results is the fact that the 
sensitivity ^t(xo) can be evaluated for trajectories within 
the attractor via consideration of the discontinuities of 
a(y). Our strategy for determining Ct(^o) is to chose the 
initial and the final separation of the trajectories to be 
the diameters Aa^o = d n . m and Axt — dn,m+t) t = 2 n — 1, 
respectively. Then, £t(xo) is obtained as 



6(zo) 



lim 

n — >oc 



(15) 



Notice that in this limit Axo — > 0, t —> oo and 
the 2 n -supercycie becomes the onset of chaos (the 2°°- 
supercycle). Then, £*(xo) can be written as 



6 (ni) 



cr n (m - 1) 



a n (m) 



, t = 2 n -l, n large, (16) 



where we have used |<7„(m)|™ ~ E[i=i \di+i,m/d^ m \ and 
rfi+i,m+2" = — di+\ >m . Notice that for the inverse pro- 
cess, starting at Axq = d n ^ m +t = —d n ,m-i and ending at 



Ax t , 



with i' = 2™ + 1 one obtains 



a n (m) 



a n (m - 1) 



f 



1, n large. (17) 



Given the known properties of <J n (m) we can readily 
extract those for £t(ra) for general non-linearity z > 1. 
Taking into account only the first 2M, M = 1,2, ... dis- 
continuities of a n (m) we have the antiperiodic step func- 
tion 



l/er n (m) 



i-M 



< m < (Z + 1)2 



i-M 



a h 12 

-ai, (2 M + l)2 n - u <m< {2 M + Z + 1)2 



when a n (m) has a discontinuity at m = l2 n ~ M . As we 
make clear below, these behaviors for the sensitivity re- 
flect the multi-region nature of the multifractal attractor 
and the memory retention of these regions in the dynam- 
ics. Thus, £t(m) = 1 (or A g (a;o) = 0) corresponds to 
trajectories that depart and arrive in the same region, 
while the power laws in Eq. 1|19[) correspond to a depart- 
ing position in one region and arrival at a different region 
and vice versa, the trajectories expand in one sense and 
contract in the other. 



IV. ORIGIN OF TSALLIS' q INDEX 

We now make explicit the mentioned link between 
the occurrence of Mori's g-phase transitions and the q- 
statistical dynamical properties at fica . Consider M = 1 , 
the simplest approximation to <r n (m), yet it captures the 
effect on £ t of the most dominant trajectories within the 
attractor. This is to assume that half of the diameters 
scales as ao = a z (as in the most crowded region of the 
attractor, x ~ 1) while the other half scales as ai = a 
(as in the most sparse region of the attractor, x ~ 0). 
With use of the identity A n = (1 + t/(2k + l))lnA/i n 2 ) 
t = 2™ — (2k + 1), £t(rn) = |ao/ai|™ can be rewritten as 
the g-exponential 



6(m)=exp 9o [A« t] , 



where 



qo = l 



In 2 



(z — 1) In a ' 



i-M 



(18) 

with I — 0, 1, 2 M — 1, and this implies that £t( TO ) = 1 
when cr„(m) is continuous at m, and 



and where 



A (fc) = 
qa {2k 



(z — 1) In a 



£t(m) 



a i 



ai+i 



or, £ t <(m) 



07 



ai+i 



(19) 



l)ln2 



k = 0,1, 



(20) 
(21) 

(22) 



The (2k + 1) term in Aq arises from the time shift 
involved in selecting different initial positions Xq ~ 1. 
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See Fig. 1 and Ref. 0. Similarly, the sensitivity for the 
trajectories in the inverse order yields 



&/(m) =exp Qo [Ag : o ) t] 



<7o 



and where A 



-2A 



(fe) 



(23) 



The factor 



q o appears because of a basic difference between 



where Q - 
of 2 in A^° 

the orbits of periods 2 n and 2°°. In the latter case, to 
reach xc ~ from xo ~ 1 at times i' = 2" + 1 the iterate 
necessarily moves into positions of the next period 2 n+ , 
and orbit contraction is twice as effective than expansion. 
Notice that the relationship between the indexes Qo = 
2 — qo for the couple of conjugate trajectories stems from 
the property exp q (x) = 1/ exp 2 _,j(— x). For z — 2 one 
obtains Qo — 1-7555 and go — 0.2445, this latter value 
agrees with that obtained in several earlier studies |j| - 

From the results for X q ^ and Aq^ wc can construct 
the two-step Lyapunov function 



A, 



(0) 



-oo < q < q , 



A(q) = i 0, g < q < 



A 



For z = 2 one has XqJ = In a/ In 2 



(0) 
Qo' 



(24) 



Qo < q < oo. 



-2A, 



(0) 
go 



1.323 and A^ = 



-2.646. See Fig. 2a. 



When the next discontinuities of importance in a(y) 
are taken into account new information on £*(a;o) is ob- 
tained in the form of additional pairs of g-exponentials 
as in Eqs. (|20() and i|23|l . In the next approximation 
for o~ n (m) there are four scaling factors, ao — ct z , oti, 
<^2, 0^3 = a, where a\ = 1/<t(1/4) (with a\ ~ 5.458 for 
z = 2) and a 2 = 1/^(3/4) (with a 2 ^ 2.195 for z = 2). 
These last values are associated to the two 'midway' re- 
gions between the most crowded and next most sparse 
regions of the attractor. We have now three values for 
the q index, qo, qi and q 2 (together with the conjugate 
values Qo = 2 — qo, Qi = 2 — qi and Q 2 = 2 — q 2 for 
the inverse trajectories). For each value of q there is a 
set of g-Lyapunov coefficients running from a maximum 
X qj .max to zero (or a minimum Xq 6 >m \ n to zero) . From the 

results for X q ^ and Xq } , j = 0, 1, 2 we can construct three 
Lyapunov functions Aj(q), — oo < q < oo, each with two 
jumps located at q = qj = 1 — In 2 / hx a j(z)/a(z) and 
q = Qj = 2 — qj. Similar results are obtained when more 
discontinuities in a n (m) are taken into account. 

Direct contact can be established now with the formal- 
ism developed by Mori and coworkers and the g-phase 
transition reported in Ref. ^4|. Each step function 
for A(q) can be integrated to obtain the spectrum </>(q) 
(A(q) = dcj)/dq) and from this its Legendre transform 
ip(X) {= <j) — {l — q)A). We illustrate this with the a n (m) 
approximated with only two scale factors and present 
specific values when 2 = 2. The free energy functions 
4>(q) and tjj(X) obtained from the two-step A(q) deter- 



mined above and shown in Fig. 2a are given by 



4° } (q - 9o), q < Qo, 



(q) 



o. 



Ag(q-Qo), q > Qo, 



qo < q < 



(25) 



and 



= I (1 - Q )A, Ag < A < 
1 (l-go)A, 0< A<A^. 



(26) 



See Fig. 2b. The constant slopes of ipW represent the 
g-phase transitions associated to trajectories linking two 
regions of the attractor, x ~ 1 and x ~ 0, and their 
values 1 — go and go — 1 correspond the index go obtained 
for the g-exponentials £t in Eqs. I|20() and 123fl . The 
slope go — 1 — —0.7555 coincides with that originally 
detected in Refs. 0|, E3- When we consider also the 
next discontinuities of importance in a(y) at y = 1/4 
, 3/4 we obtain a couple of two g-phase transitions for 
each of the three values of the g index, go, gi and q 2 . 
The constant slope values for the g-phase transitions at 
1 — go and go — 1 appear again, but now we have two other 
pairs of transitions with slope values 1 — gi and q\ — 1, 
and, 1 — q 2 and q 2 — 1, that correspond, respectively, to 
orbits that link the most crowded region of the attractor 
to the 'medium crowded' region, and to orbits that link 
this 'medium crowded' region with the 'medium sparse' 
region of the attractor. 



V. EQUALITY BETWEEN g-LYAPUNOV 
COEFFICIENTS AND RATE OF g-ENTROPY 
CHANGE 

Next, we verify the equality between the g-Lyapunov 
coefficients A^ and the g-generalized rates of entropy 

(k) 

production K qi . We follow the same procedure as in Ref. 
|6| . Consider a large number TV" of trajectories with initial 
positions uniformly distributed within a small interval of 
length Axo containing the attractor point x . A partition 
of this interval is made with N nonintersecting intervals 
of lengths £i.o (i = 1,2,..., N). For Axo sufficiently small, 
after t = 2 n — 1, or t' = 2™ + 1, iterations these interval 
lengths transform according to 



£i,0 



ai+i 



or, 



£i,0 



Oil 



ai+i 



(27) 



(Recall a n (m) has a discontinuity at m = I2 n ~ ). We 
observe that the interval ratios remain constant, that 
is, Ei.o/Axo = £i,t/Aa;t, since the entire-interval ratio 
Axt/Axo scales equally with t. Thus, the initial number 
of trajectories within each interval A/e^o/ 'Axo remains 
fixed in time, with the consequence that the original 
distribution stays uniform for all times t < T, where 
T — > oo as Ax — > 0. We can now calculate the rate 
of entropy production. This is more easily done with 
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1.323 



KT 



-1.323 



-2.646 



m - 



-A 



(a) 

v(q)- 



54 2 



-37 



21 



-1 9o =0.2445 &=1.7555 3 4 
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FIG. 2: g-phase transitions for the main discontinuity in o n (m) with index values go and Qo = 2 — go . The solid lines are the 
piece-wise continuous functions given in the text while the solid circles are the same functions calculated from the partition 
function Z(t,q) with the piece-wise as input. See text for details. 



the use of a partition of W equal-sized cells of length 
e. If we denote by Wt the number of cells that the en- 
semble occupies at time t and by Axt the total length 
of the interval these cells form, we have Wt = Axt/s 
or W t — (Ax t /Axo)(Axo/e), and in the limit e — > 0, 
Axq/e — > 1 we obtain the simple result Wt = £t- As the 
distribution is uniform, and recalling that Eq. Q for £j 
holds in all cases, the q-entropy is given by 

S q .(t)=]n qj W t = \%h, (28) 

while K ( gf=X^\ as W t=0 = 1. 

It is important to clarify the circumstances under 
which the equalities A, = K q ® are obtained as these 
could be interpreted as shortcomings of the formal- 
ism. First of all, the rate K q does not generalize the 
trajectory-based Kolmogorov-Sinai (KS) entropy /Ci that 
is involved in the well-known Pesin identity \\ = K,\ [jj - 
[Tol| . Presumably, the (/-generalized KS entropy JC q would 
be defined in the same manner as JC\ with the use of S q 
in place of Si . The rate K q is determined from values of 
S q only at two different times The relationship be- 
tween JC\ and K\ has been investigated for several chaotic 
maps [HI and it has been established that the equality 
K,\ = K\ occurs during an intermediate stage in the evo- 
lution of the entropy Si (t) , after an initial transient de- 
pendent on the initial distribution of positions and before 
an asymptotic approach to a constant saturation value. 
Here we have seemingly looked into the analogous inter- 
mediate regime in which one would expect K q = K q , as 
we explain below, however, the answer to this question 
is not analyzed in this occasion. 

We have only considered initial conditions within small 
distances outside the positions of the Feigenbaum attrac- 
tor and have not focused on the initial transient behav- 
ior referred to in the above paragraph. As for the final 
asymptotic regime mentioned above it should be kept in 
mind that the distance between trajectories, from which 
we obtain A g , always saturates because of the finiteness 
of the available phase space (the multifractal subset of 
[—1,1] that is the Feigenbaum attractor). It is widely 



known jlOj that special care needs to be taken in deter- 
mining Ai from a time series to avoid saturation due to 
folding and similar limitations occur for \ q k ^ . Separation 
of incipiently chaotic trajectories, just as separation of 
chaotic ones, undergo two different processes, stretching 
which leads to the ^-exponential regime in £ t and folding 
which keeps the orbits bounded. Therefore for t suffi- 
ciently large Eq. J7J would be no longer valid, just like 
the exponential £ t of chaotic attractors. This is the rea- 
son there is a saturation time T in our determination of 

(k) 

X q and this consequently supports our use of the rates 



VI. CONCLUDING REMARKS 

Our most striking finding is that the dynamics at the 
onset of chaos is constituted by an infinite family of 
Mori's g-phase transitions, each associated to orbits that 
have common starting and finishing positions located at 
specific regions of the attractor. Each of these transi- 
tions is related to a discontinuity in the a function of 'di- 
ameter ratios', and this in turn implies a ^-exponential 
£t and a spectrum of q-Lyapunov coefficients - equal to 
the Tsallis rate of entropy production - for each set of 
attractor regions. The transitions come in pairs with 
specific conjugate indexes q and Q = 2 — q, as these cor- 
respond to switching starting and finishing orbital po- 
sitions. Since the amplitude of the discontinuities in a 
diminishes rapidly, in practical terms there is only need 
of evaluation for the first few of them. The dominant dis- 
continuity is associated to the most crowded and sparse 
regions of the attractor and this alone provides a very 
reasonable description of the dynamics, as found in ear- 
lier studies Q - p The special values for the Tsallis 
entropic index q in £ f are equal to the special values of 
the variable q in the formalism of Mori and colleagues at 
which the g-phase transitions take place. Therefore, we 
have identified the cause or source for the entropic index 
q observed at the Feigenbaum attractor. 
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We found that the sensitivity to initial conditions at 
the onset of chaos does not have the form of a sin- 
gle g-exponential but of infinitely many interlaced q- 
cxponentials. More precisely, we found a hierarchy of 
such families of interlaced g-exponentials. An intricate 
state of affairs that befits the rich scaling features of a 
multifractal attractor. This dynamical organization is 
difficult to resolve from the consideration of a straight- 
forward time evolution, i.e. starting from an arbitrary 
position xq within the attractor and recorded at every 
time t. In this case what is observed 0] are strongly 
fluctuating quantities that persist in time with a tangled 
pattern structure that presents memory retention. On 
the other hand, if specific initial positions with known 
location within the multifractal are chosen, and subse- 
quent positions are observed only at pre-selected times, 
when the trajectories visit another region of choice, a 
well-defined g-exponential form for £ t emerges. The spe- 
cific value of q and the associated Lyapunov spectrum \ q 
can be clearly determined. For each case the value of q is 
given by the values of the trajectory scaling function a at 
one of its discontinuities, while the corresponding spec- 
trum X q reflects all starting positions in the multifrac- 
tal region where the trajectories originate. We remind 
the reader that the results presented here are indepen- 
dent from the dynamics of approach to the attractor as 
we have not considered the time evolution of initial po- 
sitions xq outside the attractor and leave this case for 
future attention. 

Interestingly, the crossover from g-statistics to BG 
statistics can be observed for control parameter values 
in the vicinity of the onset of chaos, /i > /i^, when 
the attractor consists of 2™ bands, n large. The Lya- 
punov coefficient Ai of the chaotic attractor decreases 
with Afi = fi— [Loo as Ai oc 2~ n ~ A\i K , n — ln2/ln<5(z), 
where 6 is the second Feigenbaum constant ||. The 
chaotic orbit consists of an interband periodic motion of 



period 2™ and an intraband chaotic motion. The expan- 
sion rate y~l ■ ~ In | df,, (xj)J dxj | grows as lni for t < 2" 
but as t for t » 2 n H, ]l3|. This translates as Tsallis 
dynamics with q ^ 1 for t < 2™ but BG dynamics with 
q = 1 for t > 2™. 

In summary, we have obtained further understand- 
ing about the nature of the dynamics at the onset of 
chaos in logistic maps. We exhibited links between orig- 
inal developments, such as Feigenbaum's a function and 
Mori's g-phase transitions, with more recent advances, 
such as ^exponential sensitivity to initial conditions Q 
and ^-generalized identity between Lyapunov coefficients 
and rate of entropy change 0. Chaotic orbits possess 
a time irreversible property that stems from mixing in 
phase space and loss of memory, but orbits within criti- 
cal attractors are non-mixing and have no loss of mem- 
ory. Our results apply to many other unimodal one- 
parameter families of maps. One example is the expo- 
nential function map f(y) = 1 — f exp(|y|~ z ) studied in 
Ref. as the simple (monotonic) change of variable 
\x\ z = exp(\y\~ z ) indicates. Comparable results would 
be expected to hold for the two other routes to chaos 
, intermittency and quasiperidicity, in low-dimensional 
maps as these exhibit both g-phase transitions and 
g-sensitivity to initial conditions 0|. Clearly, our 
results establish a new feasible numerical scheme to de- 
termine the family of values for the entropic index q as- 
sociated to a critical multifractal attractor (here via the 
diameter function a). Conversely, the computation of the 
sensitivity £ t at the onset of chaos offers an original means 
to evaluate Feigenbaum's trajectory scaling function a or 
its equivalent for other critical attractors. Additionally 
and remarkably, the perturbation with noise of this at- 
tractor brings out the main features of glassy dynamics 
in thermal systems poj. 
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